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Abstract — This paper addresses the stabili- 
sation of discrete-time switching linear systems 
(DTSSs) with control inputs under arbitrary 
switching, based on the existence of a common 
quadratic Lyapunov function (CQLF). The au- 
thors have begun a line of work dealing with 
control design based on the Lie-algebraic solv- 
ability property. The present paper expands 
on earlier work by deriving sufficient condi- 
tions under which the closed-loop system can 
be caused to satisfy the Lie-algebraic solvabil- 
ity property generically, i.e. for almost every 
set of system parameters, furthermore admit- 
ting straightforward and efficient numerical im- 
plementation. 

Keywords — Switching systems, Lie alge- 
bras, Common eigenvector assignment, Trans- 
verse subspaces, Genericity. 

1 INTRODUCTION 

This paper considers control design for feedback sta- 
bilisation of DTSSs under arbitrary switching regimes. 
Several results for the analysis of stability of switch- 
ing systems under arbitrary switching exist [101 [§], 
although most of these deal only with autonomous 
switching systems, i.e. switching systems without con- 
tinuous control inputs. 

The study of control design methods to achieve 
closed-loop stability under arbitrary switching for sys- 
tems with control inputs, in contrast, has been rel- 
atively scarce. Amidst the existing work, a com- 
putationally appealing approach consists in adapting 
LMI-based numerical stability tests to the synthesis 
of feedback controls that will guarantee the existence 
of a CQLF for the closed-loop switching system (e.g. 
[U [7]). Such LMI-based methods for control design, 
however, provide little information on the structure of 
the designed closed-loop system, and thus in general 
offer limited support for analysis. 

A structurally-based control design method has 
been developed by the authors [3j HI [5] that involves 
the Lie-algebraic solvability property. A well-known 



sufficient condition for stability [H] states that an au- 
tonomous switching linear system admits a CQLF 
(and is hence stable under arbitrary switching) if every 
subsystem is stable and the Lie algebra generated by 
the subsystem matrices is solvable. Solvability of a ma- 
trix Lie algebra is equivalent to the existence of a single 
similarity transformation that transforms each matrix 
into upper triangular form. The results of (3j[4jEl thus 
"activate" the aforementioned stability analysis result 
into a control design technique. 

A central contribution in [3J is an iterative design al- 
gorithm that searches for a set of stabilising feedback 
matrices that attain the target simultaneously triangu- 
larisable closed-loop structure via the application of a 
common eigenvector assignment (CEA) procedure and 
state dimension reduction at each iteration. The main 
theoretical result in [3J establishes that the proposed 
algorithm will be successful until the state dimension 
is reduced to 1 if and only if feedback matrices exist so 
that the corresponding closed-loop subsystem matrices 
are stable and simultaneously triangularisable, i.e. if 
and only if feedback matrices exist so that the closed- 
loop system satisfies the aforementioned Lie-algebraic 
stability condition. Also in [3] , a numerical implemen- 
tation for the proposed iterative design algorithm and 
CEA procedure are provided. A key structural condi- 
tion also is provided which, when satisfied, guarantees 
a directly computable solution for the CEA procedure. 
If this structural condition is not satisfied, then the 
required quantities are sought by means of an optimi- 
sation problem. 

The aforementioned Lie-algebraic stability condi- 
tion is (a) restrictive and (b) non-robust, in the sense 
that (a) it is satisfied for a very limited number of 
autonomous switching systems and (b) even if it is 
satisfied for a given system, it is almost surely not sat- 
isfied by systems with parameters arbitrarily close to 
the given one. The work in [?] then provides a ro- 
bust result by relaxing, for single input systems, the 
simultaneous triangularisation requirement to approx- 
imate (in a specific sense) simultaneous triangularisa- 
tion. The main theoretical contribution in 4 cstab- 



lishes that if a system satisfying the aforementioned 
Lie- algebraic condition exists in a suitably small neigh- 
bourhood of the given system data, then the proposed 
algorithm is guaranteed to find feedback matrices so 
that the corresponding closed-loop DTSS admits a 
CQLF even if the Lie-algebraic condition is not met 
by the given system data. 

Even if the aforementioned Lie-algebraic condition is 
restrictive for autonomous switching systems, the ex- 
istence of feedback controls causing the corresponding 
closed-loop system to satisfy such Lie-algebraic con- 
dition need not be such a restrictive condition. The 
restrictiveness of this condition is related to the key 
structural condition provided in [3]: if such structural 
condition is satisfied at every iteration of the algo- 
rithm, then the problem may be not restrictive at all 
for systems with the given dimensions. In this regard, 
the main result in [5] is the identification of the situa- 
tion that prevents the structural condition from hold- 
ing at every iteration of the algorithm. 

In the present paper, we build upon the results of 
[5] by providing sufficient conditions for the structural 
condition to hold at every iteration of the algorithm for 
almost every set of system parameters with the given 
dimensions. We thus provide sufficient conditions for 
the genericity of the property of existence of feedback 
matrices so that the closed-loop subsystem matrices 
are stable and generate a solvable Lie algebra. 

Notation. The index set {1,2, ... ,N} is denoted 
N . The kernel (null space) of a matrix or linear map 
A is denoted ker A, its image (range), imgA, and 
its spectral radius, p{A). For x £ C" xm , its trans- 
pose is denoted x' , its conjugate transpose x* and its 
Moore- Penrose generalised inverse x\ If S, T are vec- 
tor spaces, then S C T means that S is a subspace of 
T and d(«S) denotes the dimension of S. 

2 PROBLEM FORMULATION 

Consider the DTSS 

x k+ i = A i{k) x k + B t{k) u k {k) , (1) 

where x k G R n and u\ € R m * for all k, i(k) takes 
values in N_ for all k, the matrices Ai £ R ,ixn and 
Bi £ M nxm » are known for all i £ N, Bi have full 
column rank, and (Ai, Bi) is controllable for all i £ N. 
We are interested in state-feedback control design of 
the form 

u %(k) _ x i(k)Xk , (2) 
so that the resulting closed-loop system 

x k+ i = A'^ k) x k , where (3) 
Af = A i +B i K i , forzejV, (4) 

admit a CQLF and hence be stable under arbitrary 
switching. Note that at every time instant k, the con- 
trol law (|2| requires knowledge of the "active" subsys- 
tem given by i(k). 



As is well-known, ensuring that p(A° h ) < 1 for i £ 
N is necessary but not sufficient to ensure the stability 
of the DTSS ([3]) for arbitrary switching. A sufficient 
condition is given by the following result, which is a 
minor modification of Theorem 6.18]. 

Lemma 1 (Lie-algebraic-solvability stability condi- 
tion). If p(A^ L ) < 1 for i £ N_, and the Lie algebra 
generated by {A° h : i £ N_} is solvable, then 0) admits 
a common quadratic Lyapunov function and hence is 
exponentially stable. □ 

In this paper, we specifically consider stabilising 
state feedback design for the DTSS ((T|) based on the 
Lie-algebraic-solvability condition of Lemma [TJ and 
thus focus on the DTSS class defined next. 

Definition 1 (SLASF). A set Z = {(Ai £ R nxn ,B t £ 
R nxmi ) : i £ N} is said to be SLASF (Solvable Lie 
Algebra with Stability by Feedback) if there exist Ki £ 
R miXn such that A° L as in generate a solvable Lie 
algebra and satisfy p(A^ h ) < 1. □ 

In matrix terms, the fact that the Lie algebra gen- 
erated by the matrices A^ L is solvable is equivalent to 
the existence of an invertible matrix T £ C™ xn such 
that T^Af-T is upper triangular for i £ N. That 
is, each matrix A^ L is similar to an upper triangular 
matrix under a common similarity transformation T. 
Note that even if the matrices A° h have real entries, 
those of T may be complex [2] . 

3 PREVIOUS RESULTS 

Control design that causes the closed-loop system to 
be stable by satisfying the conditions of Lemma [1] can 
be performed iteratively by seeking feedback matri- 
ces that assign a common eigenvector with stable cor- 
responding eigenvalues, and reducing the state-space 
dimension by 1 at every iteration [5J [3] . This method- 
ology is given in pseudocode below as Algorithm [1] 
Algorithm [T] seeks feedback matrices Ki so that the 
closed-loop matrices given by (H]) are stable and 
simultaneously triangularisable. 

3. 1 The Algorithm 

Algorithm [1] begins by setting internal matrices equal 
to the subsystem matrices of the DTSS to be stabilised 
(Aj = Ai and B\ = Bi at the Initialisation step). At 
every iteration [I indicates iteration number, see (|5|)]. 
the algorithm executes Procedure CEA [see ©] on its 
internal system matrices (the latter matrices are A\ 
and Bf). Procedure CEA aims to compute a vector, 
v\, and corresponding feedback matrices, Ff, so that 
v\ is a feedback-assignable unit eigenvector common 
to all internal subsystems, with corresponding stable 
eigenvalues. That is, if Procedure CEA is successful, 
then v{ will satisfy \\v[\\ = 1 and (Aj+BjFf)vf = \fv{ 
for some scalars Xf satisfying |Af | < 1, for all i £ N. 
Algorithm [T] then computes internal closed-loop matri- 
ces [A i ' ch in (0)], updates internal feedback matrices 



Algorithm 1: Iterative triangularisation 
Data: A, £ R nxn , B t £ R™*™. for i E N 
Output: Ki for i E jV 
begin Initialisation 

L 4 1 = Ai, Bj = B u K? = 0, U x = I, I «- ; 
repeat 

(5) 

[«U*?}£d <- CEA({Af}f =1 ,{i?f}f =1 ) (6) 

(7) 



(8) 



(9) 

(10) 

(11) 
(12) 



if I < n then 

Construct a unitary matrix: 

Kl^l •••!<] ec-x-. 

Assign <- [«2|"-I u nj) 

<- C^iBf, 



until £ = n; 



[Kf in @] and then reduces the internal state dimen- 
sion by 1. This reduction occurs at ([9])- (fT2]) [ng is the 
internal state dimension, see ©]. Note that v\ is the 
first column of the unitary matrix , and considering 



(fTU|) then U; +1 U t+1 = I and Ug* +1 v{ = 0. Algorithm Q] 
iterates until the internal state reaches dimension 1. 
If the given system matrices form a SLASF set (re- 
call Definition [TJ , the matrices Ki computed by Algo- 
rithm [1] will be the required feedback matrices. 

If the given system matrices Ai, Bi, for i £ N_, form 
a SLASF set, then at every iteration of Algorithm [1] 
a stable feedback-assignable common eigenvector v\ 
is ensured to exist for the internal system with ma- 
trices Af, Bf, for i £ N_. Conversely, if a feedback- 
assignable common eigenvector vf exists at every iter- 
ation of Algorithm [TJ then the given system matrices 
form a SLASF set. The latter constitutes the main 
theoretical result that underpins our iterative control 
design algorithm [6J [3] . 

3. 2 The Procedure 

As expressed in the previous paragraph, the existence 
of a feedback-assignable common eigenvector with cor- 
responding stable eigenvalues is central to our devel- 
opment. This section recalls the structural condition 
introduced in [3] which, when satisfied, ensures that 
such a vector exists and allows its computation in a 
numerically efficient and straightforward way. 

We introduce some notation required to state the 
aforementioned structural condition. Define m 



rank(Bf) = d(imgS|), and factor Bf = bfrf, where 
rf : M. m * -> R m * has full row rank and bf : R m * ->• 
M™ 4 has full column rank. We adopt the convention 
that bj is an empty matrix if mf — 0. Note that 
img_Bf = imgfrf. Let A 1 be the vector with compo- 
nents Af, i £ N, i.e. 



A — [Ajl, A 2 , . . . , Ajy]', 

and build the matrix 



' X[I - A{ ■ 

Re(A e ) = : 

\ e T - A 1 



(13) 



(14) 



B l = b\kdia,g[b{,...,b i N ]., 



where blkdiag denotes block diagonal concatenation. 



Lemma 2 (Structural condition [3J[5]). Let 



N 

pe = m + ^ mf - Nng 

i=l 



(15) 



Then, 



(a) A vector that can be assigned by feedback as a com- 
mon eigenvector with corresponding eigenvalues Af 
for i £ N exists if and only if d(ker Q^(A )) > 0. 

(b) If Qi(A )w = with w ^ partitioned as 

w==[v,Ui,...,u'jj], then v ^ 0, and (16) 
(Aj + BjFf)v = \iv, for i EN, (17) 

for every Ff satisfying rfFfv = it,. For each i £ 
N one such Ff always exists and is given by Ff — 
(rf) t w i t> t . 

(c) d(ker Qe(A e )) > p£ for every choice of A 1 as in 
113\) . Consequently, if pg > 0, then a feedback- 
assignable common eigenvector exists for every 
choice of corresponding eigenvalues. □ 

Lemma [5] gives a structural condition, namely pg > 
0, for a feedback-assignable common eigenvector v to 
exist for every choice of corresponding eigenvalues Af. 
This condition is structural because the quantities in- 
volved in the computation of pg are only matrix ranks 
and dimensions. If the structural condition pg > 
is satisfied, a feedback- assignable common eigenvector 
v\ , as required at iteration I of Algorithm [TJ can be 
computed by (a) selecting its closed-loop eigenvalues 
Af corresponding to each subsystem, (b) finding a vec- 
tor w ^ with components partitioned as in (|16|) so 
that Qg(A )w — 0, i.e. so that w E kerQ^(A £ ), (c) tak- 
ing the first ng components of w, i.e. the subvector v 
in (fT6"]) . and (d) computing v[ — v/\\v\\. The feedback 
matrices that assign such eigenvector with correspond- 



= ing eigenvalues Af can be obtained as Ff = (rfyi 



Procedure CEA (Structural condition satisfied) 

Input: Af G R n ' xn \ Bf G K"* xm «, for i G N 
Output: v[, Ff for i G JV 
Factor fif = fcfrf with fof G R"^ xm * and 
mf = rank(£?f ) ; 

if p£ = ri£ + Yli=i m i ~ > then 
Select Af G M so that |Af | < 1 ; 
Find «i^0 such that Qe(A e )w = ; 
Partition w as in 02 

«i = «/IMI ; 



Ff = (rf)t Uj; wt, for i € N ; 



An implementation of Procedure CEA is thus given 
below for the case when the structural condition of 
Lemma [H is satisfied. 

Even if the DTSS matrices Ai, E>i have real en- 
tries, those of the matrices Af, Bf internal to Algo- 
rithm [Tj can be complex at some iteration £. This is 
so because the vector vf returned by Procedure CEA 
(a feedback-assignable common eigenvector) can have 
complex components even if Af, Bf have real entries, 
causing A t +1 , Bf +1 to have complex entries. However, 
when the structural condition pg > is satisfied, the 
closed-loop eigenvalues A £ can be arbitrarily selected. 
Hence, selecting real closed-loop eigenvalues will cause 
the vector v\ to have real components. In the sequel, 
we assume that real eigenvalues will be selected and 
hence all matrices internal to Algorithm [T] will have 
real entries. 

3. 3 The Structural Condition 

If the structural condition given by Lemma [3J namely 
pe > 0, holds at iteration £ of Algorithm [TJ then Pro- 
cedure CEA can easily compute a feedback-assignable 
common eigenvector and the corresponding feedback 
matrices, for every choice of corresponding closed-loop 
eigenvalues. The quantity pi depends on mf, the rank 
of Bf. At the first iteration of Algorithm [1] i.e. when 
I = 1, the internal matrices B\ = B. L have n = ni 
rows, Tfii columns, and since by assumption they have 
full column rank, then m\ — m.i. At subsequent iter- 
ations, the matrices Bf have rii = n — £ + 1 rows and 
rrii columns. According to ©-QUI) an d (H2), we have 



to* - 1 < TO- 



<TO 



(18) 

will depend 

on the feedback-assignable eigenvector v\ returned by 
Procedure CEA: 



According to ©-(HO]) and (12 



mf +1 



if v\ ^ img Bf 
1 if v { G img Bf 



(19) 



From pop , then to^ +1 = mf — 1 when mf = ng, be- 
cause necessarily in this case v\ G W 1 ' = imgBf. The 
following theorem and its corollary follow from ([T9| 
and were presented in [5]. 



Theorem 1. Consider Algorithm]]] at iteration £ and 
pn as in \15\) , with mf = rank(_B^). Then, pi+i > 
Pi — 1, with equality if and only if 



v{ G p| ivag Bf. 



(20) 



Corollary 1. Let pg > 0. Then, 

(a) p q >0 for q = £,...,£ + p e - 1. 

(b) Pi+i > if vf imgSj? for some k G N_. 

(c) Pi+i ^ if and only if pi — 1 and \20\) . 

4 MAIN RESULTS 

In this section, we derive conditions to ensure that the 
structural condition pi > will hold for £ — l,...,n. 
We will achieve this goal by looking more deeply into 
the condition (f2"0")) . In Section B. 1[ we recall a property 
of subspaces that is required for the derivation of our 
main results in Section B. 21 

4. 1 Transversality of Subspaces 

We next recall the property of transversality of sub- 
spaces (see, e.g., Chapter of [P2]). 

Definition 2 (Transverse). Two subspaces S,T of an 
ambient space X are said to be transverse when the 
dimension of their intersection is minimal, i.e. when 

d(S n T) = max{0,d(5) + d(T) - d(X)}. (21) 

Equivalently, S and 7~ are transverse when the di- 
mension of their sum is maximal. We extend this 
definition to sets of subspaces as follows. Let S = 
{S\, . . . , Sn} be a set of subspaces of an ambient space 
X . We say that S is transverse when both the intersec- 
tion of the subspaces in every subset of S has minimal 
dimension and the sum of the subspaces in every subset 
of S has maximal dimension. 

The following properties of subspaces can be 
straightforwardly established. 

Lemma 3. Let S — {Si, . . . , Sn} be a set of subspaces 
of the ambient space X , and define 

p = d(X) + ^2 d(Si) - Nd{X). 



Then, 



(a) d{Si n = d(Si) + d(Sj) - d(Si + S,). 

(b) If S is transverse, then d ^HieAr — max{0,p}. 

(c) If S is transverse and p > 0, then d(St + Sj) = 
d(X) for all i,j G N_ with i ^ j. 

(d) Let J = I U {j}, with J C N and #J = #1 + 
1. Suppose that p > and that {Si : i G /} is 
transverse. Then, {Si : i G J} is transverse if and 
only if f|ig/ Si + Sj = X. 



Proof of Lemma\3$c$. (=>) For a set K C N_, define 
p K = d(X) + E i&K d ( S i) - #Kd(X). Since p K = 
p > 0, then pi > and pj > because d(<S,) < d(X) 
for all i e N_. By Lemma [2HEJ) and since pi > and 
p, 7 > 0, then d(f] ieI Si) = pi and d(n ieJ <S;) = p,/. 
By Lemma [31Jg), we have 



d(f) S t ) = d(f) St) + d(5 3 -) - d(f) 5i + 5,) 

= p / + d(5 J -)-d(f|5j+5 i ). (22) 



is/ 



Pj 



Necessity is established by substituting the expressions 
for pi and pj into (|2"2")l and recalling that #J = #7+ 1 . 

Let K C I. We have d(A') = d(f| ieJ + < 
dCflieiir - d (' Y )- Taking if = {fc}, jointly with 
the fact that {5; : i G 7} is transverse, establishes that 
the dimension of the sum of the subspaces in every 
subset of {Si : i G 7} has maximum dimension. Also, 
we have 



d( p| S t n 5,) = d( f) 5,) + d(£,-) - d( fl S t + Sj), ^^0^ Therefore, rank^fe 1 



ieii 



ieK 



Lemma [5] below relates the dimension of the sub- 
space Sf to the controllability indices of (Af, Bf). 

Lemma 5. The dimension d(Sf) = p\ equals the num- 
ber of controllability indices equal to 1 of (Af, Bf). 

Proof. According to the standard construction for the 
controllability indices of a system (see, e.g. [T^] ), it fol- 
lows that the number of controllability indices equal to 
1 of (Af,Bf) is given by 2mf — ra.nk[/3f,Af/3f], where 
(3f is any matrix satisfying img (3f — img Bf . Since 
Sf C Bf, write Bf = Bf © Sf and let a = d(Bf). 
Then, p\ = m\ — a. Let {b\, . . . ,b a { be a basis 
for Bf, {b a+ i, . . . ,b m e} be a basis for Sf, and j3f — 
[6i,...,6 m |]. By Lemma 14©, A% £ Bf for k = 

1, . . . , a and Afbk G Bf for k = a + 1, . . . , mf. There- 
fore, rank^Af^f] < mf + a. If rank[/3f , Af(3f] < 

i 

mf + a, then J2T=i c fii + J2t=i dkAfbk = for some 
scalars Cj and dk, where not all the dk are zero. Then, 
A iY,k=i d kbh G Bf and AfJ2k=i d k b k £ Sf, a con- 

a and 
□ 



2mf 



d(A-) 



which, jointly with the fact that {Si : i G 7} is trans- 
verse, establishes that the dimension of the intersec- 
tion of the subspaces in every subset of {Si : i G J} 
has minimum dimension. □ 

As is well known [T3] , the property of transversality 
is generic, i.e. it is satisfied for almost every set S 
composed of a finite number of subspaces of X selected 
"randomly" among all subspaces of X. 

4. 2 Genericity of the SLASF Property 

As previously mentioned, we will look more deeply into 
the condition (|20p. We define the following 



Our main result is given below as Theorem [5] We 
will provide comments and explanations after its proof. 
The proof of Theorem [5] requires an additional result, 
given as Lemma |B1 

Theorem 2. Let {Sf : i G N} be transverse, and 
qi > 0. Then, pg > for I = 1, . . . , n. 

Lemma 6. Consider Algorithm]]] at iteration £. Sup- 
pose that (Af,Bf) is controllable and Af^vf = Xfvf 
with v\ 7^ and scalar Xf. Then, <S i +1 D U^ +l Sf, 
(A^ 1 , Bf +1 ) is controllable, and 



B< 



J+ 1 



(23) 



> 



1 if v{e St, 

otherwise. 



(27) 



iEN 



According to (|23|. then (j20|) can be rewritten as v[ G 
B e . Recall that v\ is a feedback- assignable common 
eigenvector. The following result is straightforward. 

Lemma 4. Let Sf be the set of vectors v G Bf for 
which there exist Ff and X with \X\ < 1 so that 

(Af + BjFf)v = Xv. (24) 

(a) The set Sf is a subspace. 

(b) v G Sf if and only if v G Bf and Afv G Bf. 

By definition, Sf is the set of feedback-assignable 
eigenvectors for the subsystem (Af,Bf) that are con- 
tained in Bf. Consequently, vf G Bf if and only if 
v{ G Sf . In the sequel, we will employ the following. 



Proof. Let {tj : j = 1, 



if} be a basis for Bf and 



let K%j, for j — l,...,mf be the controllability in- 
dices of (Af,Bf). By ([7]) and the feedback invariance 
of controllability indices, nf • also are the controlla- 
bility indices of the pair (A 1 ^ , Bf) . Since (Af,Bf) is 
controllable, then (A^ clj ,Bf) also is controllable, and 

D = {(Af Ch )Hj :, 1 /»;:/■• ^ - 1} 

is a basis for M. ne . Write v\ with respect to the basis 
D: v{ — J2j k c j,k{Ai CL ) k tj, where not all the Cj t k are 
zero. Combining the latter with Af ch v{ = Xfv[ yields 



■jAA ) tj = 2^ A * c i' fc( A ) tj 



(28) 



3,k 



s l 



d(Sf), 

n^f. 



q e = ne + ^2pj- Nnt, (25) 
ieN_ 

p«±d(S*). 



(26) 



From (|2"5|) . it follows that Cj.k ^ for at least one pair 
of indices (j,k) such that jfe = k| • — 1, or otherwise 
the vectors in D would be linearly dependent, a con- 
tradiction. Let R — maxj{«;f ■ : Cj t k with k = 
K ij ~ 1}) an d let I be such that c-^ R -\ ^ 0. From 



the basis D, construct another basis, D, by replac- 
ing the basis vector (A e i ' CL ) K, ~ 1 t I by v[. Note that 
Span{J/; +1 £ : t G D} = and U* t+1 v[ = 

(recall Section l3~Tjl. By ©-flTTJ and the fact that 
Af 0h v{ = \jv{, then U* l+1 Af L = ^+ x t/; +1 . Hence 

^; +1 (A- ,m )%- = (^ +1 )*t7+i*i«««i 

{(4 +1 ) fc C/; +1 ^:j = l,...,mf; 

fc = 0,...,4 i -l,0',fc)^(i,K-l)} (29) 

is a basis for R™ €_1 [recall that, by ([5]), ni + i = ni~ 1]. 
From |Q2J), it follows that = f/; +1 Sf. We have 

that a basis for is E = {U£ +1 tj : j = 1, . . . , mf } 

if k > 1 or E = {U; +1 tj : j = l,...,mf;j ^ 1} if 
k — 1 . The condition k = 1 hence happens if and 
only if v\ 6 6>,f . The preceding derivations show that 
the controllability indices of (A* +1 ,Bf +1 ) are given 
by nfj 1 = for j — l,...,mf with j ^ I and 
Atj't 1 = k| r — 1 whenever k — n\ l > 1. From the lat- 
ter expressions, and recalling Lemma [SJ ([2"T]) and the 
controllability of (A^ +1 ,B e t +1 ) are established. Note 
that pl +1 = p\ + 1 whenever k = 2, since then 
= k — 1 = 1 and hence {Af +1 ,B% +1 ) has one 
controllability index equal to one more than (Af, Bf). 
The fact that Sf +1 D U£ +1 Sf follows from the latter 
consideration and the basis E. □ 

Proof of Theorem^ First, we prove that, if true, the 
following conditions 

{Sf : i G N_} transverse, qg > 0, 
(Al,Bf) controllable, 

imply that pg > 0. Since Sf C Bf, then p\ < mf and 
qi < pg. From controllability of (A\ , Bf) and Lemma[5l 
then p\ — m\ if and only if m\ — ng. Hence, if qg = pg, 
then pg — ng > 0. Otherwise, < qi < pg. 

Next, we establish the validity of (|30p for £ = 
1, . . . , n. Note that (|30[) hold at I = 1 by assumption. 
Next, suppose that (|30|) hold at some 1 < £ < n— 1. By 
the argument in the previous paragraph, then pg > 0, 
which ensures the existence and computation of v\ =^ 
such that Al' ch v{ = \\v\ with scalar \\ for all i € N. 
Hence, ( J 4f +1 ,Sf +1 ) is controllable by Lemmai Also 
by Lemma H we have S l t +1 D U£ +1 Sf for all i G JV. 
Since qg > 0, by Lemma [3 Tfb|) we have that p = qg, 
and from Lemma we have d(«S| + Sj) = ng for all 
i,j £ N with i ^ j. It follows that d(Sf +1 + > 
d(C/; +1 (5f + 5j)) = n t - 1 = n £+i for all *, j G jV 
with i / j. The latter fact establishes that the sum 
of the sets in every subset of {S- +1 :i£N} has max- 
imal dimension and also that {<Sf +1 , is trans- 
verse for all i, j G N_ with i ^ j. Let T be a subset of 
{Sf +1 : i G JV}. We proceed by induction on the num- 
ber of subspaces in T. We have already established 
that T is transverse if #T = 2. Suppose next that T 
is transverse whenever #T = 2, . . . , a, with a < N—l. 



Let T = : i G /}, with / C JV and #/ = a, and 

consider R — T U so that #R = a + 1. By 

Lemma O and properties of maps and subspaces, we 
have 

^U; +1 (f]S f i+ S e 3 ). (31) 
iei 

By fl30|) and since I C N, then {Sf : i G /} is 
transverse. By Lemma G3|d]), then d(f] ieI S- + Sj) = 
ng. Combining the latter equality with (1311) . then 
d(n ie /<Sf +1 +Sj +1 ) =ng-l = n l+1 . By Lemma [30 
then R is transverse. We have thus established that 
our induction hypothesis is valid for a + 1 and we 
conclude that {S^ +1 : i G N} is transverse. By 
Lemma [3T|bt . then p e+1 — ma,x{0,qg + i}. From ([23)1 
and (1271) . it follows that the minimum value for qg+\ 
is qi — 1, and this happens only if = pf — 1 for all 
i £ N. However, if g<> = 0, then p^ +1 > p\ for at least 
one i G N_ because, since p e = d(S e ) = qg = 0, then 
v\ £ S . Consequently qg+i > and hence we have 
established §U§ for I = 1, . . . , n. □ 

Theorem [5] gives a condition, namely {S} : i G N} 
transverse and q\ > 0, under which the structural 
condition of Lemma [TJjnj) is satisfied at every itera- 
tion of Algorithm [T] The quantity q\ depends on the 
dimensions of the subspaces S} for i G N, which, 
by Lemma [5] equals the number of controllability in- 
dices equal to 1 of (A\ — A^B] = Bi). A con- 
sequence of Corollary 5.4 of [TJ] is that the latter 
number is generically (i.e. for almost every matri- 
ces Ai and Bi) nonzero when im > n/2 and if the 
latter holds, generically equal to m, — (n mod rrii). 
If > n/2, arbitrary choices for the entries of Ai 
and Bi yield arbitrary Sf, although generically of di- 
mension rrii — (n mod rrii). Therefore, if the system 
dimensions n and m, for i G N_ are such that q\ > 0, 
then {S} : i G N} will be transverse generically in 
the space of parameters of the matrices Ai and J3,, for 
i G N_. For example, a DTSS with two subsystems 
(N = 2), order n — 6, subsystem 1 having 4 inputs 
(mi = 4) and subsystem 2 having 5 inputs (rrii = 5) 
will generically satisfy the hypotheses of Theorem [2j 
since p\ = mi — (n mod mi) = 2, p\ — 4 and hence 
qi = 6 + (2+4) — 2-6 = 0. From the preceding analysis, 
it follows that the hypotheses of Theorem [5] can hold 
only for DTSSs where each subsystem has "a lot of" 
inputs (m, > n/2, q\ > 0). 

5 CONCLUSIONS 

We have addressed feedback stabilisation of discrete- 
time switched linear systems with control inputs. The 
control strategy employed is to seek feedback matrices 
so that the closed-loop subsystem matrices are stable 
and generate a solvable Lie algebra. The problem of 
feedback stabilisation by means of the latter strategy 



is known to not always have a solution. In this context, 
we have derived conditions under which this problem 
is ensured to have a solution for most possible sets of 
parameters. These conditions hold only for DTSSs of 
specific system dimensions, where each subsystem has 
a considerable number of inputs, as compared with the 
system dimension. 
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